Mon. Not. R. Astron. Soc. 0Q0,[Tin3](2Q09) Printed 19 February 2009 (MN KTeX style file v2.2) 



A downturn in intergalactic C IV as redshift 6 is approached 
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ABSTRACT 

We present the results of the largest survey to date for intergalactic metals at redshifts z > 5, 
using near-IR spectra of nine QSOs with emission redshifts Zom > 5.7. We detect three strong 
ClV doublets at z^hs — 5.7-5.8, two low ionisation systems at Zabs > 5, and numerous 
Mg 11 absorbers at Zabs = 2.5-2.8. We find, for the first time, a change in the comoving mass 
density of CiV ions as we look back to redshifts z > 5. At a mean {z) = 5.76, we deduce 
^^c IV — (4.4 ±2.6) X 10^^ which implies a drop by a factor of ~ 3.5 compared to the value 
at z < 4.7, after accounting for the differing sensitivities of different surveys. The observed 
number of C IV doublets is also lower by a similar factor, compared to expectations for a 
non-evolving column density distribution of absorbers. 

These results point to a rapid build-up of intergalactic C IV over a period of only 
^ 300 Myr; such a build-up could reflect the accumulation of metals associated with the 
rising levels of star formation activity from z ^ 9 indicated by galaxy counts, and/or an 
increasing degree of ionisation of the intergalactic medium (IGM), following the overlap of 
ionisation fronts from star-forming regions. If the value of VIqiv we derive is typical of the 
IGM at large, it would imply a metallicity Zjgm ^ lO^"* Zq. The early-type stars responsible 
for synthesising these metals would have emitted only about one Lyman continuum photon 
per baryon prior to z — 5.8; such a background is insufficient to keep the IGM ionised and 
we speculate on possible factors which could make up the required shortfall. 
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1 INTRODUCTION 

It is quite remarkable that elements heavier than helium exist in the 
intergalactic medium (IGM) only ~ 1 Gyr after the Big Bang, as 
indicated by the presence of metal absorption lines associated with 
the Lya forest of high redshift QSOs (Songaila 2001; Becker et 
al. 2006; Ryan-Weber, Pettini & Madau 2006; Simcoe 2006). For 
these elements to be observed there needs to be significant star for- 
mation prior to redshift 2 = 6, and the nucleosynthetic products of 
such early stars need to be expelled beyond the gravitational influ- 
ence of their host galaxies. The production of early metals is closely 
linked to the emission of the Lyman continuum (LyC) photons 
responsible for ionizing the IGM and ending the so-called "dark 
ages" (Gnedin & Ostriker 1997; Fan, Carilli, & Keating 2006), a 
process which could have started as early as 2; = 20 and continued 
until completion at 2 ~ 6 (Dunkley et al. 2009). 

Clearly, some galaxies were already actively forming stars at 

2 > 6. The current record for a spectroscopically confirmed galaxy 
is held by the Lya emitter at z — 6.96 discovered by lye et al. 
(2006), but several candidates selected via the redshifted Lyman 
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break have been found at 2 > 7 (Stark et al. 2007; Richard et al. 
2008; Bouwens et al. 2008 and references therein). Furthermore, 
analyses of the broad-band spectra of galaxies at redshifts 2 = 2 — 6 
have in several cases revealed objects with 'mature' stellar popula- 
tions, suggesting that star formation in these galaxies was underway 
well before 2 = 6 (Shapley et al. 2001; Shapley et al. 2005; Yan et 
al. 2006; Eyles et al. 2007). 

Whether the level of star formation so far accounted for at 
2 > 6 is in fact sufficient to maintain reionization is far less certain, 
because it depends on a number of essentially unknown parameters, 
especially the escape fraction of Lyman continuum photons from 
the sites of star formation, the degree of clumping of the sources 
and of the IGM, and the slope of the faint end of the galaxy lumi- 
nosity function at these early times (Madau, Haardt, & Rees 1999). 
Using the luminosity function of Bouwens et al. (2006) at 2 ~ 6 
and assuming an escape fraction of 20%, Bolton & Haehnelt (2007) 
concluded that there may be just enough LyC photons to keep the 
IGM ionised, suggesting that cosmic reionization was a 'photon 
starved' process (see also Gnedin 2008). QSOs are unlikely to help 
matters, since their contribution to the metagalactic photoionization 
rate at these redshifts is estimated to be at most 15% (Srbinovsky & 
Wyithe 2007). Bolton & Haehnelt (2007) proposed that, in order to 
complete reionization by 2 ~ 6, either the escape fraction of LyC 
photons must have been higher at earlier times, or the ionizing spec- 
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trum must have been harder (or both). AUematively, the faint end 
slope of the luminosity function may have been steeper at 2 > 6 
compared to the luminosity function determined by Bouwens et al. 
(2006), as indeed proposed by Stark et al. (2007). 

The cosmic budget of elements produced via stellar nucle- 
osynthesis (loosely referred to as metals) provides an entirely dif- 
ferent, and complementary, census of the star-formation activity 
prior to a given epoch. Indeed, at redshifts 2 = 2 — 3 there ap- 
pears to be an approximate agreement (within a factor of about two) 
between the observed metal budget — summing together the contri- 
butions of stars, IGM and diffuse gas in galaxies — and expectations 
based on the integral of the cosmic star formation rate density over 
the first 2-3 Gyr of the Universe history (Pettini 2006; Bouche et al. 
2007). 

At higher redshifts, we rely increasingly on the Lya forest 
and gamma-ray burst afterglows to trace metals in the IGM and 
in galaxies, since at 2 > 5 the spectra of even the most luminous 
galaxies are too faint to be studied in detail with current instru- 
mentation. The most commonly encountered metal lines associated 
with the Lya forest are the C IV AA1548.2041, 1550.7812 doublet 
(Cowie et al. 1995; Ellison et al. 2000). Photoionisation modelling 
(Schaye et al. 2003; Simcoe, Sargent, & Ranch 2004) suggests that 
at redshifts 2 = 2 — 3 the carbon abundance is a function of the 
gas density, as expected (Cen & Ostriker 1999), and that the mean 
metallicity of the IGM is Ziqm ~ 1/1000 Zq. 

Interestingly, the comoving mass density of triply ionised car- 
bon, obtained by integrating the distribution of column densities 
N{C IV), appears to be approximately constant from 2 ~ 5 to 
2 ~ 1.5: expressed as a fraction of the critical density, Q,cn — 
(1 - 4) X 10"** (Songaila 2001; Pettini et al. 20o£l) over a pe- 
riod of ~ 3 Gyr which saw the peak of the cosmic star-formation 
and massive black hole activity (e.g. Reddy et al. 2008 and refer- 
ences therein). This surprising result has led to conflicting interpre- 
tations as to the source of these intergalactic metals, with some au- 
thors (e.g. Madau, Ferrara, & Rees 2001; Porciani & Madau 2005) 
proposing an origin in low-mass galaxies at 2 ^ 6, when the pol- 
lution of large volumes of the IGM would have been easier, while 
others (e.g Adelberger 2005) place them much closer to the massive 
galaxies responsible for most of the metal production at redshifts 
2 = 1.5-5. While both processes presumably contribute, their rel- 
ative importance has yet to be established with certainty (Songaila 

2006) . 

The most extensive theoretical work in this area has been pub- 
lished by B. D. Oppenheimer and collaborators. In a series of 
papers (Oppenheimer & Dave 2006, 2008; Dave & Oppenheimer 

2007) these authors analysed the output of cosmological simula- 
tions of galaxy formation which include feedback in the form of 
momentum-driven winds. They found that the apparent lack of evo- 
lution of Q.C IV from 2 = 5 to 1.5 can be explained as the result of 
two countervailing effects: an overall increase of the cosmic abun- 
dance of carbon, reflecting the on-going pace of star formation, and 
an accompanying reduction in the fraction of carbon which is triply 
ionised. 

Whether this is the true picture or not, it is clear that de- 
termining Q.C IV at 2 > 5 is the next observational priority; in 
any case, metal lines are our only remaining probe of the IGM 
at these high redshifts where the Lya forest itself becomes effec- 
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lively opaque. Such an extension of metal-line surveys in QSO 
spectra necessitates observations at moderate spectral resolution 
and signal-to-noise ratio at near-infrared (near-IR) wavelengths, a 
regime which has only recently begun to be exploited for such pur- 
poses (e.g. Kobayashi et al. 2002; Nissen et al. 2007). In a pilot 
study (Ryan- Weber, Pettini, & Madau 2006), we demonstrated that 
QSO absorption line spectroscopy could indeed be performed in 
the near-IR to the levels required to detect intergalactic CiV ab- 
sorption. From observations of two QSO at redshifts 2om = 6.28 
and 5.99 we discovered two strong C IV doublets at absorption red- 
shifts 2abs = 5.7238 and 5.8290. These two absorbers, which were 
confirmed by an independent study by Simcoe (2006), would — if 
typical — imply that a high concentration of C IV was already in 
place in the IGM at 2 ~ 6, but clearly much better statistics than 
those afforded by only two QSO sight-lines are necessary to assess 
the true level of Q.c iv- 

Given the success of our pilot study, in the last two years we 
have conducted an observational campaign on the Keck and ESO- 
VLT telescopes aimed at securing near-IR spectra of as many as 
possible of the 19 known QSOs with 2cm > 5.7 and infrared mag- 
nitude J < 19.2. The results of this programme are presented 
here. Specifically, we have observed 13 QSOs and obtained useful 
data for ten of them. This sample represents an increase by a fac- 
tor of seven in the pathlength probed, from the absorption distance 
AX = 3.6 covered by Ryan- Weber et al. (2006) to AX = 25.10 
With this extended data set, we find the first evidence for a decrease 
in Q.C IV from 2 ~ 4.7 to 2 ~ 5.8, in broad agreement with the pre- 
dictions of the cosmological simulations by Oppenheimer, Dave, & 
Finlator (2009). 

This paper is organised as follows. Section 2 describes the 
observations, followed by brief comments on each QSO in Sec- 
tion 3. In Section 4 we deduce the value of Q,cn implied from 
our three detections of C IV absorbers within the redshift range 
5.2 ^ 2abs ^ 6.2. In Sections 5 and 6, we consider the complete- 
ness of the survey, and discuss the evidence for a drop in Q.c n as 
we look back to redshifts 2 > 4.7. The implications of these results 
for early star formation and reionisation are discussed in Section 7. 
We finally summarise our main conclusions in Section 8. 

Throughout this paper we use a '737' cosmology with Hq = 
70 km s~^ Mpc"^, f^M = 0.3andf2A = 0.7. Much of the previous 
work (including our own) with which we compare the present re- 
sults was based on an Einstein-de Sitter cosmology; in all cases, we 
have converted the quantities of interest to the cosmology adopted 
here. 



2 OBSERVATIONS AND DATA REDUCTION 

A trawl of the relevant literature reveals that 19 QSOs are currently 
known with emission redshift 2cm > 5.7 and J magnitude brighter 
than J ~ 19.2. This magnitude is at the faint limit for record- 
ing medium-resolution near-IR spectra of S/N ^ 5 with integra- 
tion times of a few hours on 8-10 m class telescopes using avail- 
able instrumentation. The redshift cut is required to give a useful 
pathlength over which to search for C IV absorption at redshifts 
-Zabs > 5, this being the limit of previous surveys (Songaila 2001; 
Pettini et al. 2003). 

^ The absoiption distance X is defined so that dX/dz = Hq{1 + 
z)'^ /H(z), where H is the Hubble parameter. Populations of absorbers 
with constant physical cross sections and comoving number densities main- 
tain constant line densities dN/ dX as they passively evolve with redshift. 
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Table 1. Details of the Observations 



QSO 


(mag) 




Telescope/ 
Instrument 


Wavelength 
Range (A) 


Resolution Integration 
(kms"-*^) Time (s) 






ULAS J()20332.38+001229.2 


19.1 


5.706<= 


Keck II/NIRSPEC 


10 219-10 588 


185 


6 300 


10-13 


5.601-5.639'' 


SDSS JOS 1 827.40+ 1 7225 1 . 8 


18.5 


6.00 


VLTl/ISAAC 


9 877-10 790° 


53 


32400'" 


8-13 


5.380-5.930'' 


SDSS J083643. 85+005453.3 


17.9 


5.8108 


VLTl/ISAAC 


10100-10 550 


53 


16 200 


8 


5.524-5.752 


SDSS J084035.09+562419.9 


19.0 


5.774<= 


Keck II/NIRSPEC 


9 548-11 140 


185 


11700 


10-22 


5.167-5.706 


SDSS J103027.01+052455.0h 


18.9 


6.309' 


VLTl/ISAAC 


9 880-10 780= 


53 


82 800j 


5-7 


5.382-5.951« 


SDSS Jl 13717.73+354956.9 


18.4 


5.962<^ 


Keck II/NIRSPEC 


9 542-11 172 


185 


9 000 


20^0 


5.163-5.893 


SDSS Jl 14816.64+525150.2 


18.1 


6.421'^ 


Keck II/NIRSPEC 


9 572-11 159 


185 


18 000 


16-30 


5.183-6.196 


SDSS J130608.26+035626.3'^ 


18.8 


6.016' 


VLTl/ISAAC 


9 914-10 332 


53 


7 200 


5 


5.404-5.662 


SDSSJ160253.98+422824.9 


18.5 


6.051<= 


Keck II/NIRSPEC 


9 553-11 143 


185 


9900 


17^8 


5.170-5.981 


SDSS J205406.49 -000514.8 


19.2 


6.062 


Keck II/NIRSPEC 


9681-10927 


185 


5 400 


8-12 


5.253-5.991 


''Signal-to-noise ratio over wavelength range sampled. 
''Redshift range covered for C IV absorption. 

'^When available, we quote values of Zcm measured from the C IV A1549.062 broad emission line in our NIRSPEC spectra. As is well known, 
the QSO systemic redshift, ^sys, can be higher than the value measured from high ionisation lines such as C IV, by up to several thousand km 



(e.g. Richards et al. 2002). 

''BAL QSO; not included in the survey statistics. 

°With a small gap between two wavelength settings. 

'Two wavelength settings, with exposure times of 14 400 and 1 8 000 s respectively. 
EKurk et al. (2007). 

''Observations reported in Ryan- Weber et al. (2006). 
'Jiang et al. (2007). 

Jin three, partially overlapping, wavelength settings, with exposure times of 39 600, 7 200 and 36 000 s respectively, 
''^em measured from the Si IV A1396.752 broad emission line in our NIRSPEC spectrum. 



Over the four-year period between 2004 and 2008 we ob- 
served 13 of these 19 possible targets with NIRSPEC on the 
Keck II telescope (McLean et al. 1998) and with ISAAC on 
VLT-UTl (Moorwood et al. 1998). For ten QSOs we secured 
spectra with S/N 5; relevant details are collected in Ta- 
ble [T] In the other three cases, the object was not pursued 
after initial test exposures, either because the spectrum dis- 
played broad absorption lines (SDSS J104433.04-012502.2 at 
z^^ = 5.74 and SDSS J104845.05+4637 18.3, Zom = 6.20), 
or because the QSO turned out to be fainter than anticipated 
(SDSS J141 11 1.29+121737.4, Zcm = 5.93). The values of emis- 
sion redshift Zcm listed in Table[T]are those measured from the dis- 
covery spectra, unless otherwise noted. 

2.1 NIRSPEC Observations 

With NIRSPEC we used the NIRSPEC 1+thin blocker filter com- 
bination to record the wavelength region ~ 9540-11 170 A on the 
ALADDIN-3 InSb detector which has 1024 x 1024 27^im pixels. 
With the 42 x 0.57arcsec entrance slit, the spectrograph delivers 
a resolution FWHM = 6.4 A (185 km s"^ at 10355 A), sampled 
with 2.1 pixels. Typically, we exposed on an object for 900 s before 
reading out the detector, moving the target by a few arcseconds 
along the slit, and starting the subsequent exposure. Total integra- 
tion times ranged between 5400 and 18 000 s (see column (7) of Ta- 
ble[B. All of the QSO observed were sufficiently bright at near-IR 
wavelengths to be visible on the slit, thereby facilitating guiding. 

The NIRSPEC two-dimensional (2-D) spectra were processed 
with a suite of purpose-designed IDL routines kindly provided by 
Dr. G. Becker. As described in more detail by Becker, Ranch, & 
Sargent (2009), this software does a very good job at subtracting 
the dark current, removing detector defects, flat-fielding, and par- 
ticularly removing the sky background (achieved by employing the 



optimal sky subtraction techniques of Kelson 2003) which can be 
the limiting factor to the S/N achievable at near-IR wavelengths. 
Optimal extraction (Home 1986) is also used to trace the QSO sig- 
nal in the processed 2-D images and output a vacuum-heliocentric 
wavelength calibrated one-dimensional (1-D) source spectrum and 
associated error. The wavelength calibration uses the OH sky emis- 
sion lines superposed on the QSO spectrum, with the wavelengths 
tabulated in the atlas by Rousselot et al. (2000). 

As we have many (^ 6) individual 900 s spectra (and asso- 
ciated error) of each of our QSOs, we combined them using a 
weighted mean algorithm that operates on a pixel-by-pixel basis 
and rejects outliers beyond a specified number of standard devia- 
tions from the mean (usually ^ 3ct). The spectra were finally nor- 
malised by dividing by a spline fit to the underlying QSO contin- 
uum. The rms deviation of the data from this continuum (away from 
obvious spectral features) provides an empirical estimate of the fi- 
nal S/N achieved0 We found this empirical estimate of the S/N to 
be lower, by factors between 2.0 and 1.1, than the pixel-to-pixel 
error on the weighted mean calculated when co-adding the indi- 
vidual exposures, as described above. This difference presumably 
reflects the residual importance of systematic (rather than random) 
errors; in any case, the values listed in the penultimate column of 
Table[T]are those measured from the residuals about the continuum 
fit. Since the S/N varies along each spectrum, we give in the Table 
the maximum and minimum values which apply to the wavelength 
region searched for absorption lines (listed for each QSO in column 
(5) of TablelB. 



Strictly speaking, this is a lower limit to the real S/N of the spectra, since 
it assumes a 'perfect' fit to the QSO continuum. However, in practice this is 
the value to be used in assessing the significance of any absorption features 
present in the spectrum. 
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Figure 1. NIRSPEC spectraof ULAS J020332.38+001229.2, SDSS J084035.09+562419.9 and SDSS J113717.73+354956.9. In each panel, the top histogram 
shows the normalised QSO spectrum, while the lower histogram is the corresponding Icr error spectrum. Absorption features are indicated by vertical tick 
marks and labelled with their identification. The two CiV absorption doublets seen in the spectrum of ULAS J020332. 38+001229. 2 (top panel) are close to 
the QSO emission redshift and are probably intrinsic (see section lXTl . Parameters of the C IV doublet identified towards SDSS J113717. 73+354956. 9 (bottom 
panel) are listed in Tables|2]and[3] 



2.2 ISAAC Observations 

The ISAAC spectra were obtained with a mixture of service- 
mode and visitor-mode observations We used the instrument in its 
medium resolution mode with a 0.6 arcsec wide slit, which results 
in a spectral resolution FWHM = 1.8 A (53 km s~^), sampled 
with four pixels of the 1024 x 1024 Hawaii Rockwell array. With 
this set-up, the wavelength range spanned is AA ~ 460 A; we nor- 
mally chose a grating setting which placed the C IV emission line 
at the red limit of the range covered. In some cases, multiple wave- 
length settings were used, to increase the pathlength for absorption 
covered (see TablelTJ. The observations were conducted in the con- 



ventional 'beam-switching' mode, in a series of 4 x 900 s exposures 
at two locations on the detector, typically separated by 10 arcsec. 
After each cycle of four integrations, the object was moved along 
the slit to a new pair of positions. The instrument was aligned on 
the sky at an angle chosen to include on the slit a bright star which 
provides a useful pointing and seeing reference. 

The reduction of the ISAAC data used mainly IRAF tasks fol- 
lowing the steps already described by Ryan- Weber et al. (2006). 
Briefly, we started from the ESO pipeline products which are 2-D 
images that have undergone a substantial amount of pre-processing, 
including flat-fielding, background subtraction, and wavelength 
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Figure 2. NIRSPEC spectra of SDSS Jl 14816.64+525150.2, SDSS J160253.98+422824.9 and SDSS 1205406.49-000514.8. In each panel, the top histogram 
shows the normalised QSO spectrum, while the lower histogram is the corresponding Icr error spectrum. Absorption features are indicated by vertical tick 
marks and labelled with their identification. 



calibration. Given the weakness of the sources at the moderately 
high dispersion of the ISAAC spectra, we found it most effective 
to combine all the processed frames of a given QSO (at a given 
wavelength setting) in 2-D, after appropriate shifts to bring them 
into alignment. This co-addition also allowed us to reject deviant 
pixels. 

One-dimensional spectra were extracted from the co-added 2- 
D frames, with a second pass at subtracting any residual sky fea- 
tures; the data were then rebinned to two wavelength bins per res- 
olution element in order to improve the S/N while still maintaining 
adequate spectral sampling. The final steps, including normalisa- 
tion to the QSO continuum and estimate of the rms noise, were the 
same as described for the NIRSPEC spectra. As can be appreci- 



ated from inspection of Table [T] the values of S/N achieved with 
ISAAC tend to be lower than those of the NIRSPEC spectra, de- 
spite the generally longer integration times. However, the signifi- 
cantly higher spectral resolution provided by ISAAC, by a factor of 
~ 3.5, largely compensates for the lower S/N, so that the two sets 
of data reach comparable sensitivities in terms of the minimum ab- 
sorption line equivalent widths detectable. 



3 NOTES ON INDIVIDUAL QSOS 

The normalised spectra of the ten QSOs in Table[T]are reproduced 
in Figures [T}0 together with the corresponding error spectra. Ab- 
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Figure 3. ISAAC spectra of SDSS J081827.40+172251.8 (two wavelength settings) and SDSS J083643.85+005453.3. In each panel, the top histogram shows 
the normalised QSO spectrum, while the lower histogram is the corresponding Icr error spectrum. Absorption features are indicated by vertical tick marks and 
labelled with their identification. There is a tentative detection of an intervening C IV absorber towards SDSS J081827. 40+172251. 8 at z^bs = 5.7899; this 
possible system was omitted from the sample of definite detections used to estimate He iv ■ 



sorption lines were identified from a visual inspection of the spec- 
tra; they are labelled in Figures [T}@] and tabulated in Table |2] We 
now briefly comment on each QSO spectrum in turn. 



3.1 ULASJ020332.38+001229.2 

This QSO was identified by Venemans et al. (2007) from the first 
data release of the UKIRT Infrared Deep Sky Survey (UKIDSS), 
and independently by Jiang et al. (2008) from the Sloan Digital 
Sky Survey (SDSS). Both surveys reported similar emission red- 
shifts: Zcm ~ 5.86 (Venemans et al.) and Zem = 5.854 (Jiang 
et al.). From our NIRSPEC spectrum, however, we measure the 



lower value Zcm = 5.706 from a broad emission feature which 
we presume to be ClvA1549.062; this value is in much better 
agreement with the recent reassessment of the optical and infrared 
spectrum of this source by Mortlock et al. (2008) who conclude 
that 2em ~ 5.72. Mortlock et al. point out that this is a broad ab- 
sorption line (BAL) QSO which led the discoverers to misidentify 
the N V A1240.1 emission line as Lya: the latter is completely ab- 
sorbed by the blueshifted N V absorption trough. 

As this is a BAL QSO, we do not include it in the survey for 
C IV absorption (as we are primarily interested in C IV in the in- 
tergalactic medium, and wish to avoid contamination of our sam- 
ple by metal-enriched gas which may be ejected from the QSO). 
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Figure 4. ISAAC spectra of SDSS J103027.01+052455.0 (two wavelength settings) and SDSS J130608.26+035626.3. In each panel, the top histogram shows 
the normalised QSO spectrum, while the lower histogram is the corresponding la error spectrum. Absorption features are indicated by vertical tick marks and 
labelled with their identification. Parameters of the two C IV doublets identified towards SDSS J103027.01+052455.0 are listed in Tables [2] and[3] 



Nevertheless, it is of interest to note that we do detect two narrow 
Civ doublets at redshifts Zabs = 5.7541 and jZabs = 5.7961 (see 
Figure [TJ; their equivalent widths are given in Table [2] Relative to 
Zcm = 5.706, these two absorbers have velocities v = +2150 and 
+4025 km s~^ and are thus most likely associated with the QSO 
environment. Positive velocities relative to z^m by up to ~ 2000- 
3000 km s^^ have been encountered before among high ionisation 
absorbers at lower redshifts (e.g. Fox, Petitjean, & Bergeron 2008; 
Wild et al. 2008); the values we deduce here are rather extreme and 
may in part reflect the uncertainty in the systemic redshift Zgys. 



3.2 SDSSJ081827.40+172251.8 



This Zcm ~ 6.00 QSO was observed with ISAAC at two 
wavelength settings which, together, cover the wavelength re- 
gion 9877-10 790 A with only a very small gap between them 
(see top two panels of Figure [Sj. We find two possible ab- 
sorption systems: a C IV doublet at ^abs ~ 5.7899 and a 
Mg II AA2796. 3543, 2803.5315 doublet at ^abs = 2.8341; their 
reality can only be confirmed with higher S/N data. 
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Table 2. Absorption Lines Identified 



QSO 






Ion 




Wo (A) 


Comments 


ULAS J()20332.38+001229.2 


5.706 


5.7540 


Civ 


1548.2041 


0.38 


Intrinsic 






5.7543 


Civ 


1550.7812 


0.25 


Intrinsic 






5.7959 


Civ 


1548.2041 


0.62 


Intrinsic 






5.7962 


Civ 


1550.7812 


0.56 


Intrinsic 


SDSSJ081827.40+172251.8 


6.00 


5.7897 


Civ 


1548.2041 


0.12 


Tentative Identification 






5.7911 


Civ 


1550.7812 


0.06 


Tentative Identification 






2.8348 


Mgll 


2796.3543 


0.18 


Tentative Identification 






2.8334 


Mgll 


2803.5315 


0.22 


Tentative Identification 


SDSS J084035.09+562419.9 


5.774 


5.5940 


Sill 


1526.70698 


0.52 








5.5938 


Fell 


1608.45085 


0.34 








5.5943 


AllI 


1670.7886 


0.44 




SDSS J103027.01+052455.0 


6.28 


5.7238 


Civ 


1548.2041 


0.65 








5.7241 


Civ 


1550.7812 


0.41 








5.8299: 


Civ 


1548.2041 




^abs 3nd Wo uncertain 






5.8299: 


Civ 


1550.7812 




•^abs '^O uncertain 


SDSS Jl 13717.73+354956.9 


5.962 


5.0127 


AllI 


1670.7886 


0.70 








5.0120 


Fell 


1608.45085 


0.65 








2.6260 


Mgll 


2796.3543 


0.52 








2.6261 


Mgll 


2803.5315 


0.55 








5.7380 


Civ 


1548.2041 


0.33 








5.7387 


Civ 


1550.7812 


0.20 




SDSS Jl 14816.64+525 150.2 


6.421 


2.5077 


Mgll 


2796.3543 


0.56 


Tentative Identification 






2.5084 


Mgll 


2803.5315 


0.29 


Tentative Identification 






2.7727 


Mgll 


2796.3543 


0.77 








2.7728 


Mgll 


2803.5315 


0.60 








3.5545 


Fell 


2344.2139 


0.34 








3.5544 


Fell 


2382.7652 


0.96 




SDSS J130608.26+035626.3 


5.99 


2.5309 


Mgll 


2803.5315 


2.50 


A2796.3543 outside spectral range'^ 


SDSS J160253.98+422824.9 


6.051 


2.7138 


Mgll 


2796.3543 


0.44 








2.7135 


Mgll 


2803.5315 


0.38 




SDSS J205406.49-0005 14.8 


6.062 


2.5973 


Mgll 


2796.3543 


2.05 








2.5995 


Mgll 


2803.5315 


2.29 





''From the compilation by Morton (2003). 

''This value of equivalent width is in the observed frame (since the feature is unidentified). 
'^Both members of the doublet are present in the GNIRS spectrum published by Simcoe (2006). 



3.3 SDSSJ083643.85+005453.3 

No absorption lines were identified in the 450 A stretch of the near- 
IR spectrum of this QSO recorded with ISAAC (see bottom panel 
of Figure |3j. 



which may be associated with neutral gas in galaxies; however, it is 
hard to envisage how this can be accomplished, given that the asso- 
ciated Lyman series lines are essentially indistinguishable against 
the backdrop of the nearly-black Lya forest. 



3.4 SDSSJ084035.09+562419.9 

Our NIRSPEC spectrum of this QSO includes the broad C IV emis- 
sion line at redshift Zc-m = 5.774, to be compared with the SDSS 
value 2om = 5.85 ± 0.02 reported by Fan et al. (2006b). We find 
no C IV systems, but we do detect three absorption lines from a low 
ionisation system at Zabs = 5.5940 (middle panel of Figure [T] and 
Table I2J. It would be of great interest to determine the metallic- 
ity (and other physical parameters) of such high redshift absorbers 



3.5 SDSSJ103027.01+052455.0 

This was the first QSO to be observed as part of this programme 
and a portion of its spectrum was published by Ryan-Weber et 
al. (2006). In Figure |4] we show the full ISAAC spectrum which 
extends from 9880 A to 10 780 A, recorded with three, partially 
overlapping, grating settings. We find two strong C IV doublets, at 
Zabs = 5.7238 and 5.8288. The latter was considered as marginal 
by Ryan- Weber et al. (2006), but independent observations by Sim- 
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Table 3. Parameters of identified C iv absorption systems 



QSO 



Instrament 



^abs 



log AfCClvVcm^^) f)(kms^i) Comments 



SDSSJ103027.01+052455.0 ISAAC 5.7238 ± 0.0001 14.39 ± 0.05 61 ± 6 

SDSSJ103027.01+052455.0 ISAAC 5.8288 ± 0.0002 14.45 ±0.13 26 ± 9 

SDSSJl 13717.73+354956.9 NIRSPEC 5.7383 ± 0.0002 14.20 ±0.15 30 ± 12 



SDSS JOS 1 827.40+ 1 7225 1 . : 



ISAAC 



5.7899 ± 0.0003 13.57 ±0.16 



39 ±25 



Tentative identification 



coe (2006) with the Gemini Near-Infrared Spectrograph (GNIRS) 
confirm its reality. Even so, it is difficult with the quality of our 
data to measure the equivalent widths of the two C IV lines in this 
absorber, as they appear to be blended with other broader features. 



3.6 SDSSJ113717.73+354956.9 

This is one of our best NIRSPEC spectra (see FigurelTJ. The emis- 
sion redshift we determine from C IV A1549.062, Zcm ~ 5.962, 
is lower than the value Zcm = 6.01 ± 0.02 reported by Fan et al. 
(2006b) from the discovery SDSS spectrum. We identify absorp- 
tion lines from three separate systems in the 1630 A portion of the 
spectrum recorded with NIRSPEC: a high redshift, low ionisation 
system at Zabs = 5.0124, a Mg II absorber at ^abs ~ 2.6261, and a 
Civ doublet at Zabs ~ 5.7383. The equivalent widths of all these 
lines are given in Table |2] 



3.7 SDSSJ114816.64+525150.2 

This is the highest redshift QSO in the present sample. We mea- 
sure Zcm = 6.421 from the Si IV A1396.752 broad emission line, 
in good agreement with the value Zam ~ 6.42 published by Fan et 
al. (2006c). Our NIRSPEC spectrum (top panel of Figure|2j shows 
a definite Mg II doublet at Zabs = 2.7728, and a possible one at 
2abs = 2.5080 which appears to be part of a broad feature of un- 
certain identity. We also cover two Fe II lines, AA2344, 2382, in 
a strong low ionisation system at Zabs = 3.5544; other Fe II and 
Mg II lines in this system can be recognised in the spectrum pub- 
lished by Barth et al. (2003)0 



3.8 SDSSJ130608.26+035626.3 

The relatively short integration (two hours) of this object with 
ISAAC resulted on a relatively noisy spectrum (S/N ~ 5) repro- 
duced in the bottom panel of Figure |4] and already reported by 
Ryan- Weber et al. (2006). No CiV doublets are detected. The 
strong absorption line present in the spectrum is the longer wave- 
length member of the Mg II AA2796.3543, 2803.5315 doublet at 
Zahs = 2.5309; both lines can be discerned in the GNIRS spec- 
tra of this object obtained independently by Simcoe (2006) and by 
Jiang et al. (2007). 



^ We are grateful to the referee, Xiaohui Fan, for pointing out the existence 
of this absorption system. 



3.9 SDSSJ160253.98+422824.9 

Our measured Zam = 6.051 is in reasonable agreement with Zam ~ 
6.07 from the SDSS (Fan et al. 2006c). We detect a definite Mg II 
system at Zabs = 2.7137, as well as a puzzling broad absorption 
feature (see middle panel of Figure |2j. If this were a C IV 'mini- 
BAL', its ejection velocity relative to Zem would be 7150km s~^. 



3.10 SDSSJ205406.49-000514.8 

Even a relatively short exposure (5400 s) with NIRSPEC is suf- 
ficient to reveal the very strong Mg II doublet present at ^abs ~ 
2.5984 (bottom panel of Figure|2j. This Mg II absorber, and the one 
at 2abs = 2.5309 in SDSS J130608.26+035626.3, are close to the 
threshold for classification as "Ultra-Strong Systems" in the SDSS 
survey by Nestor, Turnshek, & Rao (2005). It is interesting to note 
that these authors find a steep redshift evolution in the frequency of 
such systems, which become progressively rarer at 2; 1. 



4 SURVEY FOR INTERGALACTIC CIV AT Z > 5 
4.1 Redshift Range Sampled 

In the last column of Table [T] we have listed for each QSO the 
redshift intervals over which we searched for C IV absorption. The 
values of Zmin, which range from 5.16 to 5.52, correspond to wave- 
lengths close to the values of Amin in column (5) for the detection 
of A1548.2041. The values of Zmax, which range from 5.66 to 6.20, 
correspond either to wavelengths close to the values of Amax in col- 
umn (5) for the detection of A1550.7812, or to Zcm ~ 3000 km s"\ 
whichever is the lower [the QSO emission redshifts are listed in col- 
umn (3)]. We exclude 'proximate systems' within 3000km s^^ of 
the emission redshift because they may be associated with the QSO 
environment, rather than with the more generally distributed IGM; 
however, none were found in the present survey except for the two 
systems in ULAS J020332.38+00 1229.2 which is not included in 
the analysis anyway as it is a BAL QSO (see section lSTt . 

From these values of Zmin and Zmax we calculate the corre- 
sponding absorption distances, given by: 

xiz) = -^{[nM{i + zf + n^]'/^~i} (I) 

which is valid for Qm + f^A = 1, where Qm and Qa are re- 
spectively the matter and vacuum density parameters today. With 
f^M ~ 0.3 and Qa = 0.7, we find that our survey covers a total 
absorption distance AX — 25.1. 
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4.2 C IV Column Densities 

Summarising the results of section[3] we find three definite and one 
possible C IV doublets over this absorption distance. We deduced 
values of column density log A'^(C IV)/cm~^ for these absorbers by 
fitting the observed absorption lines with theoretical Voigt profiles 
generated by the VPFIT packagelf] For each absorber, VPFIT de- 
termines the most likely values of redshift Zabs, Doppler width 6 
(km s~^), and column density log A'^(ClV)/cm^^, by minimizing 
the difference between the observed line profiles and theoretical 
ones convolved with the instrumental point spread function. Vac- 
uum rest wavelengths and /-values of the C IV transitions are from 
the compilation by Morton (2003). For an illustrative use of this 
line-fitting software, see Rix et al. (2007). 

Table[3]lists the values returned by VPFIT, together with their 
Icr errors, for the four C IV absorbers considered here. Strictly 
speaking, these values of A'^(C IV) are lower limits to the true col- 
umn densities because at the relative coarse spectral resolution of 
our data the absorption lines we see are likely to be the unresolved 
superpositions of several components (e.g. Ellison et al. 2000), 
some of which may be saturated. This concerns applies particu- 
larly to the zabs = 5.7383 doublet in SDSS J113717.73-I-354956.9 
for which the b value on which VPFIT converged (b = 30 km s^^, 
to be interpreted as an 'equivalent' velocity dispersion parame- 
ter describing the ensemble of components) is much smaller than 
the NIRSPEC instrumental broadening fcinstr = llOkms"^ — in 
other words, these doublet lines are totally unresolved. On the other 
hand, the doublet ratio W(1548)/Vy (1550) = 1.65 is consistent 
with only a moderate degree of line saturation. In the other three 
cases, the equivalent b-values are comparable to, or larger than, 
= 32km s~^ of ISAAC (see Table [3}, indicating that the 
C IV absorption lines are partially resolved. For the remainder of 
our analysis, we shall assume that the column densities of CiV 
have not been underestimated by large factors, although this as- 
sumption will need to be verified in future with higher resolution 
observations. 



4.3 Mass density of C iv 

It has become customary to express the mass density of an ion (in 
this case C IV) as a fraction of the critical density today: 



IV = 



Hp mc r 

C Pcrit 



Nf(N)dN, 



(2) 



where Hq = 100 ft km s Mpc is the Hubble constant, 

fyic IV 

is the mass of a C IV ion, c is the speed of light, pcrit = 1-89 x 
10"'^^ g cm^'^, and f{N) is the number of CiV absorbers per 
unit column density per unit absorption distance X{z). In the ab- 
sence of sufficient statistics to recover the column density distribu- 
tion function f{N), as is the case here, the integral in eq. ^ can 
be approximated by a sum: 



Nf{N)dN = 



E.iV.(Civ) 



AX 



with an associated fractional variance: 



[E.^^civ)] 



(3) 



(4) 
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Figure 5. Cosmological mass density of C IV as a function of redshift. The 
blue squares show the measurements by Songaila (2001), the red triangle is 
from Pettini et al. (2003), and the green triangle is the value deduced here. 
All values plotted have been reduced to the '737' cosmology adopted in the 
present work. EiTor bars are Icr. While this plot shows the actual values of 
Sic IV measured, they are not strictly comparable because each of the three 
surveys had a different sensitivity limit. This issue is discussed in detail in 
the text (section|5). 



as proposed by Storrie-Lombardi, McMahon, & Irwin (1996). With 
h = 0.70, eqs. Q and l|3) then lead to: 



nc<v = 1.63 X 10" 



AX 



(5) 



Summing up the values of A^i(ClV) in Tableland dividing by 
AX = 25.1, we find: 



f2c,v = (4.4 ±2.6) X 10" 



(6) 



VPFIT is available from |http://www.ast.cam.ac.ukr| rfc/vpfit.html 



over the redshift interval Zabs ~ 5.2-6.2 at a mean (^abs) = 5.76. 
We have not included the possible system at Zabs = 5.7899 in 
the summation; its inclusion would increase Qcn by only 5%. For 
comparison, Pettini et al. (2003) measured ilciv ~ (2.2 ± 0.8) x 
10^** at (^abs) ~ 4.69 (in the same cosmology as adopted here); 
these two values, together with the lower redshift determinations 
by Songaila (2001), are plotted in Figure[5] 

Thus, at first sight, Qc iv seems to have dropped by about a 
factor of five in only ~ 300 Myr, as we move back in time from 
z ~ 4.7 to 2; ~ 5.8. However, before we can be confident of 
such a rapid build-up of C IV at these early epochs in the Universe 
history (these redshifts correspond to only ~ 1 Gyr after the Big 
Bang — see top axis of Figure |5j, we have to look carefully at the 
different completeness limits of the surveys whose results are col- 
lected in Figure[5] In particular, the sensitivities of absorption line 
searches at near-IR wavelengths, such as the one reported here, are 
still lower than those achieved at optical wavelengths with compa- 
rable observational efforts and we suspect that this difference must 
be contributing to some extent to the high redshift drop evident in 
Figure |5] 



5 COMPLETENESS LIMIT OF THE SURVEY 

In order to quantify this effect, we conducted a series of Monte 
Carlo simulations, as follows. We used VPFIT to generate theoreti- 
cal C IV doublets for ranges of values of b and 7V(C IV); these fake 
absorption lines were then inserted at random redshifts into the real 
spectra, after convolution with the appropriate instrumental profile 
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Figure 6. Sensitivity function to C IV doublets in our survey, deduced from 
1440 Monte Carlo realisations. See text (section|5) for further details. 



(depending on whether the spectrum to be 'contaminated' was a 
NIRSPEC or an ISAAC one) and with random noise reflecting the 
measured S/N of the spectrum at that wavelength. We then tested 
whether we could recover these C IV doublets by visual inspection 
of the spectra, in the same way as the real absorbers were identi- 
fied. So as to test the sensitivity of our survey as comprehensively 
as practical, we chose three sets of values for the velocity disper- 
sion parameter, h — 20, 40 and 60 km s^^, reflecting the range 
of fe-values of the real C IV doublets (see Table [3). For the col- 
umn density we considered twelve values, in 0. 1 dex steps between 
log A'^(C iv)/cm^^ — 13.4 and 14.5. For each combination of b 
and log7V(C iv), we created 40 Monte Carlo realisations, so that 
the full series of completeness tests involved 1440 fake C IV dou- 
blets distributed within our parameter space. 

We found that the fraction of C IV doublets recovered, g{N), 
did not depend sensitively on b nor Zahs but, as expected, was a 
function of column density — see Figure|6] While we seem to be ~ 
90% complete for column densities log Ai'(C iv)/cm~^ J5 14.2, 
our sensitivity drops fairly rapidly below this value: we can only 
detect about half of the C IV doublets with half this column density 
and, if N(C IV) drops by a further factor of two, we miss some 85% 
of the absorbers. For comparison, the estimate of fie iv by Pettini et 
al. (2003) refers to C IV absorbers with log iV(C iv) /cm"^ > 13.0 
(and includes a 22% correction for the fall in completeness of their 
Keck-ESI spectra from ~ 100% at log7V(C iv)/cm"^ ^ 13.3 to 
45% at logAf(Clv)/cm"^ = 13.0). Songaila's (2001) estimates 
of ilc IV , on the other hand, refer to the yet wider column density in- 
terval 12.0 log A'^(C iv)/cm~'^ ^ 15.0 which she could access 
thanks to the high spectral resolution and S/N of her Keck-HIRES 
spectra. From her Figure 1, it appears that her number counts be- 
come progressively less complete below log Af(C iv)/cm~^ ~ 
13.0, although no correction was applied to account for this. 

From the above discussion it is evident that, in order to assess 
the significance of the suggested drop in Qc iv beyond z ~ 4.7, we 
need to consider the results of these different surveys over the same 
range of C IV column densities. Referring to Figure |5] we see that 
we are > 50% complete for log A'^(C iv)/cm~^ > 13.8, so that 
we can compare with the surveys by Songaila (2001) and Pettini 
et al. (2003) in this column density regime. Specifically, we have 
adjusted the values flc iv reported by those authors as follows: 



C,3.8 Nf{N)dN 
C"o" iV/(iV)diV 



Figure 7. Cosmological mass density of C IV as a function of redshift, 
taking into account only CiV absorbers with column densities 13.8 ^ 
log Af(C iv)/cm~^ ^ 15.0, to which all three surveys considered here 
are sensitive. The symbols have the same meaning as in Figure m For the 
highest redshift data point in the sample by Songaila (2001), shown by the 
pale blue square, the correction of eq. (7) probably underestimates the value 
ofC'c,,. 



for the values published by Songaila (2001), and 

Nf{N)dN 



NfiN)dN 



(8) 



for the value published by Pettini et al. (2003), assuming 

f{N) = 10'"■''^iV^3'•^ (9) 

where N13 is A'^(ClV) in units of 10^^ cm~^, as determined by 
Songaila (200 ifl 

The lower limits of integration in the denominators of eqs. ^ 
and ([Sj reflect the sensitivity limits of the corresponding surveys. 
Since the survey by Songaila (2001) is only partially complete 
for Civ absorbers with 12.0 < log iV(C iv)/cm"^ 13.0, 
eq. Q may overestimate the correction required — this works in the 
sense of diluting the evidence for a drop in Qcn at higher red- 
shifts. The effect is most pronounced for the highest redshift bin 
in Songaila's sample, because her data is least sensitive to column 
densities log 7V(C iv)/cm^^ < 13.0 at Zabs ^ 4.5 (see Figure 2 
of Songaila 2001); accordingly we have indicated her highest red- 
shift data point with a different colour in Figure|7] 

For the NIRSPEC+ISAAC sample presented here, we have 
modified eq. Q to include our sensitivity function g{N) shown in 
Figure |6] 



ficiv = 1.63 X 10^ 



E,A^i(Civ) X <7(iV.) 



AX 



(10) 



(7) 



All the errors have been scaled accordingly, so as to maintain the 
same percentage error. 

Values of il'c n are plotted in Figure [7] Evidently, we still see 
a drop in the mass density of C IV as we look back beyond 2 ~ 4.7, 
albeit at a reduced significance compared to Figure|5] Considering 
only absorbers with log A''(C iv) /cm^^ ^ 13.8, we find that Q'q „ 
decreases by factor of ~ 3.5, as we move from ficiv = {^-^ ± 
0.6)xlO~*at(2) = 4.6910^0 IV = (5±3) x 10"^ at {2) =5.76 
(1(7 errors). The drop is significant at the 4ct level. 



Again, adjusted to the present cosmology. 
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6 ABSORPTION LINE STATISTICS 

We can examine the evidence for a decrease in intergalactic C IV at 
redshifts z > 4.7 in an alternative way from considering simply the 
sum of their column densities. It is instructive to ask how the num- 
ber of detected C IV doublets compares with the number we may 
have expected under the assumption of an invariant column density 
distribution, f{N). Again, we assume f{N) = lO^^^'^^Affg^ * 
(Songaila 2001), which implies /(> iV) = 2.iN^^-^ per unit 
absorption distance. 

Our survey is not uniform in either S/N nor resolution, and 
both factors determine the minimum value of N{C IV) detectable. 
We therefore calculated the expected value of /(> A*') on a spec- 
trum by spectrum basis, using as threshold the 5(t detection limit 
for the equivalent width of an absorption line: 

Wo(5cr) = X — (11) 



l + z 



1 

sTn 



where 5\ is the resolution in A. Values of resolution and S/N for 
each spectrum are listed in Table [T] With Wo(5cr) is associated a 
minimum column density: 



iVciv(5cr) = 1.13 X 10^" X 



A2/ 



(12) 



where A and Wo are both in A and / is the oscillator strength of the 
transition (in this case C IVA1550.7812, the weaker member of the 
doublet). 

In Figure [8] we plot the cumulative number of absorbers ex- 
pected (for the minimum column density to which each spectrum 
is sensitive) as a function of the cumulative pathlength. Over our 
total survey pathlength, AX — 25.1, we expect to detect lOlg 
Civ absorbers. We instead find three absorbers, indicating that 
the distribution of CiV column densities which applies at red- 
shift 1.5 ^ z ^ 4.5 does not continue unchanged to 5.2 5C 
z ^ 6.2. The likelihood of detecting three C IV systems when 
ten are expected is less than 1%. It is interesting to note that the 
absorption systems we do detect were found in spectra with col- 
umn density sensitivities in the middle of the range of our sur- 
vey: iVciv(5cr) = 9.1 X 10" cm"^ and 1.2 x lO" cm"^ for 
SDSSJl 13717.73-1-354956.9 and SDSS J103027.01-I-052455.0 re- 
spectively. The absorber towards SDSS Jl 13717.73-1-354956.9 rep- 
resents a 6.6(T detection. Thus, it seem unlikely that we are missing 
a considerable number of absorption systems above our 5(t sensi- 
tivity limit. 



7 DISCUSSION 

To summarise, by extending our previous searches for intergalactic 
C IV to near-IR wavelengths, we have found that both the number 
of such absorbers per unit absorption distance and the comoving 
mass of triply ionised carbon are approximately three times lower 
at 2 ~ 5.8 than at z < 4.7. Although our results refer to the 
stronger C IV systems in the column density distribution function, 
with log A^(C iv)/cm~'^ ^ 13.8, the evolution appears to extend 
at least to absorbers with half this column density (Becker et al. 
2009). 

Before considering the implications of these results, we reflect 
briefly on possible caveats. As already mentioned, we may have un- 
derestimated the column densities of C IV through line saturation; 
we do not expect this effect to be so severe as to cause (spuriously) 
the fall in fie iv evident in Figure[7] but only higher resolution ob- 
servations will quantify the corrections required (if any). It must 
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Figure 8. Cumulative number of C IV absorption systems expected from our 
survey (if there is no evolution in the column density distribution function 
from redshift z ~ 3 to 6). The cumulative number is plotted as a function 
of the survey pathlength, AX, which increases as each spectrum is added 
to data set. We expect to find 10 C IV absorption systems in our total path- 
length of AX = 25.1. The star symbol denotes the three absorbers we 
detect in our entire survey in the redshift range 5.2 < z < 6.2. 



also be borne in mind that the distribution of C IV at z > 5.2 seems 
to be extremely patchy: of the nine sight-lines sampled, seven show 
no absorption, while two absorbers are found in the direction of one 
QSO. Thus, our statistics are very prone to what in galaxy counts 
is referred to as 'cosmic variance'. This can explain why the pilot 
observations by Ryan- Weber et al. (2006) suggested a higher value 
of f^c iv-the first stage of our survey fortuitously found one-third of 
the NiiC iv) of the whole sample presented here, even though 
it covered only ~ 15% of the total AX. Conversely, Becker et al. 
(2009) detected no absorbers in their higher resolution NIRSPEC 
spectra of four QSOs (three in common with our sample) over less 
than half of the AX of our survey. In future, better statistics will 
undoubtedly improve the accuracy of determinations of Slciv at 
> 5; for the moment we shall proceed on the assumption that 
our estimate of the error on this quantity, as shown in Figure |7] is 
realistic. 

The downturn in S7civ we have discovered as we approach 
2: ~ 6 is clear evidence of important changes taking place in the 
Universe at, or before, this epoch. In the following discussion we 
shall find it more convenient to 'turn the time arrow around', and 
to think of the evolution as time progresses from higher to lower 
redshifts. The rise in Slciv we see in Figure |7] could result either 
from an increase in the amount of intergalactic carbon, or from 
changes in the ionisation balance favouring the fraction of carbon 
which is triply ionised, or both. Indeed, there are indications from 
previous work that both effects may be at play. 

By integrating the galaxy luminosity function down to 1/5 
of the fiducial luminosity at 2 = 3 (that is, considering galaxies 
with L ^ 0.2 L*^3), Bouwens et al. (2008) found that the cosmic 
star-formation rate density increased by a factor of 4 over the 
400 Myr period from 2 = 9 to 2 = 6. The growth in intergalactic 
C IV we see at 2 < 6 may well be a reflection of this enhanced star- 
formation activity, as the products of stellar nucleosynthesis accu- 
mulate in the IGM, presumably transported by galactic outflows. 
Note that this is carbon produced by massive stars, with evolution- 
ary time-scales of less than 10 Myr (Pettini et al. 2008b; Pettini 
2008); it is therefore quite conceivable that over ~ 400 Myr the 
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products of Type II supemovae could have travelled considerable 
distances from their production sites. 

On the other hand, Fan et al. (2006c) have argued for a marked 
increase in the transmission of the Lya forest at z ^ 6, which has 
been interpreted as marking the overlap of the hydrogen ionisation 
fronts from star-forming galaxies, as the IGM became transparent 
to Lyman continuum photons (e.g. Gnedin & Fan 2006). Opinions 
are still divided, however, as to the correct interpretation of the 
Lyof forest evolution at these high redshifts (e.g. Becker, Ranch, 
& Sargent 2007 and references therein). Coupled with the rise in 
the comoving density of bright QSOs (Fan 2006), a higher IGM 
transmission would quickly ionise a larger fraction of carbon into 
Civ. 

To make progress beyond these qualitative remarks, detailed 
modelling is necessary. As mentioned in the Introduction, B. Op- 
penheimer and collaborators have devoted considerable attention 
to this problem and shown, in particular, how momentum-driven 
winds from star-forming galaxies can reproduce the behaviour of 
f2civ at redshifts z; < 5. It is therefore of interest to consider 
their model predictions for earlier epochs. Their latest models (Op- 
penheimer, Dave, & Finlator 2009) do indeed entertain a marked 
growth in Q.c iv from z = 8 to 5, driven by both effects mentioned 
above. In their scenario, the overall density of carbon increases by 
a factor of ~ 2 between z = 5.8 and 4.7. The accompanying rise 
in f^c IV depends on the spectral shape of the radiation field which 
dominates the ionisation equilibrium. Oppenheimer et al. (2009) 
consider two possibilities: the metagalactic background in the for- 
mulation by Haardt & Madau (1996) with recent updates, and a 
more local field dominated by early-type stars in the nearest galaxy 
to the C IV absorbers, in recognition of the fact that cosmic reioni- 
sation was likely to be very patchy at these early epochs. The for- 
mer produces the larger change in Q.c iv (for the reasons explained 
above), which grows by a factor of ~ 5 between z = 5.8 and 4.7; 
the latter results in a more modest rise, by a factor of ~ 2 — 3. 

Given our warnings about the effects of 'cosmic variance' on 
the current observational determinations of Sic iv, it would be quite 
premature to compare in detail these theoretical predictions with 
the data in Figure |7] and, for example, draw conclusions on the 
shape of the ionising background at these early epochs. It is encour- 
aging, nevertheless, that models and data concur in a general trend 
towards a build-up of C IV over the relatively short cosmic period 
between 2: ~ 6 and 5. Two ingredients of the simulations by Op- 
penheimer et al. (2009) need to be clarified, however. The first, and 
more important one, is that in their scenario the galaxies responsi- 
ble for most of the IGM enrichment are of low mass and with lumi- 
nosities well below the current observational limit L ^ 0.2 L't^-i,. 
We know very little, if anything at all, about the galaxies at the faint 
end of the luminosity function at 2 > 5, and this situation is un- 
likely to be remedied until the launch of the James Webb Space 
Telescope in the next decade. A second concern is the resonant 
opacity of the IGM at the frequencies of the Hell Lyman series 
which has been neglected until recently, and which may affect the 
distribution of carbon among its highly ionised states (Madau & 
Haardt 2009). 

7.1 Reionization 

The average metallicity of the IGM at a given epoch reflects the 
accumulation of the heavy elements synthesised by previous gen- 
erations of stars and expelled from their host galaxies. The mas- 
sive stars which explode as Type II supemovae and seed the IGM 
with metals are also the sources of nearly all of the Lyman contin- 



uum photons produced by a burst of star formation. Consequently, 
it is relatively straightforward to link a given IGM metallicity with 
a minimum number of LyC photons which must have been pro- 
duced up until that time. The close correspondence between the 
sources metals and photons makes the conversion from one to the 
other largely independent of the details of the stellar initial mass 
function (IMF). For example, Madau & Shull (1996) calculate that 
changing the index of a power-law formulation of the IMF, such as 
Salpeter's (1955), from a = 2 to 3, results in only a 10% difference 
in the conversion between metallicity and number of LyC photons 
emitted. 

When in the past this line of reasoning has been applied to 
measures of the IGM metallicity at 2 ~ 3 (e.g. Madau & Shull 
1996; Miralda-Escude & Rees 1997), it has been concluded that 
the star formation activity prior to that epoch probably produced 
sufficient photons to reionise the Universe. With our new estimate 
of ficiv, much closer in time to what may have been the end of 
the process of reionisation, we are able to revisit this question more 
critically. 

We begin by calculating the IGM metallicity by mass, Ziqm, 
implied by our measure of Sic iv : 



f^Civ Ctot 



(13) 



rib C IV Ac 

where f7b = 0.0224//i^ (Pettini et al. 2008a) is the contribution of 
baryons to the critical density, C IV/Ctot is the fraction of carbon 
which is triply ionised, and Ac is the mass fraction of metals in 
carbon. For comparison, in the Sun, = 0.0122 and Ac ~ 
0.178 (Asplund, Grevesse, & Sauval 2005); here we assume that 
the same value of Ac applies at high redshifts because: (a) the ratio 
of carbon to oxygen (the major contributor to Zq ) is approximately 
solar at low metallicities (Pettini et al. 2008b; Fabbian et al. 2009); 
and (b) any departures are in any case likely to result in only a 
small correction compared to other uncertainties in the following 
discussion. We also adopt the conservative limit Civ/Ctot ^ 
0.5, as done by Songaila (2001); this is the maximum fractional 
abundance reached by C IV under the most favourable ionisation 
balance conditions. With these parameters, our measured 



flciv = (4.4 ±2.6) X 10" 



(14) 



implies 

ZiGM = (1.1 ± 0.6) X 10"'' = (9 ± 5) X lO"'^ Zq. (15) 

In recent years, much theoretical attention has been directed to esti- 
mating the likely properties of metal-poor early-type stars. Schaerer 
(2002) pointed out that stellar luminosities at far-UV wavelengths 
are higher at lower metallicities, leading to a greater output of LyC 
photons associated with stellar nucleosynthesis. Specifically, the 
energy emitted in hydrogen ionising photons (per baryon) is related 
to the average cosmic metallicity by: 



Ez 



r)mpC {Z) MeV 



(16) 



where mpC = 938 MeV is the rest mass of the proton and 77 is 
the stellar efficiency conversion factor. Schaerer (2002) estimated 
?7 = 0.014 for stars with Z = 1/50Zq, a factor of ~ 4 higher 
than at solar metallicities (see also Venkatesan & Truran 2003). 
For an average energy of 21 eV per LyC photon (Schaerer 2002), 
our determination of {Z\cm) ~ 1.1 x lO^'' then implies that less 
than one (~ 0.7) LyC photon per baryon was emitted by early-type 
stars prior to 2 = 5.8. 

This background is clearly insufficient to keep the Universe 
ionised. When we take into account that: (a) an unknown fraction. 
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/osc < 1, of the emitted LyC photons will escape from the re- 
gions of star formation into intergalactic space, and (b) more than 
one intergalactic LyC photon per baryon is likely to be required 
to keep the clumpy IGM ionised against radiative recombinations 
(e.g. Madau, Haardt, & Rees 1999; Bolton & Haehnelt 2007), we 
reach the conclusion that our shortfall is by at least one order of 
magnitude. 

There are a number of possible solutions to this problem. 
The most obvious one is to recall that eq. l ll5t gives the lower 
limit to the IGM metallicity, for two reasons. First, our account- 
ing of fJciv only includes high column density absorbers, with 
A'^(C Iv) ^ 10^** cm~^. A steepening of the column density distri- 
bution, f{N), at z > 5.2 could hide significant amounts of metals 
in lower column density systems. However, the recent null result 
by Becker et al. (2009) runs counter to this explanation. 

Second, if C IV/Ctot 0.5, some of the tension would be 
relieved. For example, in the models considered by Oppenheimer 
et al. (2009) most of the C is doubly-ionised al z — 5-6, and 
C IV/Ctot ~ 0.1, a factor of five lower than the upper limit as- 
sumed in eq. dlSl l. On the other hand, there are limits to the ex- 
tent to which the shortfall can be explained by ionisation effects 
alone. If the IGM is not highly ionised at these redshifts, so that 
Civ/Ctot <C 0.5, we would also expect a non-negligible frac- 
tion of C to be singly ionised. Unlike C III, which is virtually im- 
possible to observe at these redshifts because its resonance lines 
are either too weak or buried in the Lya forest, C II has a well- 
placed strong transition at rest wavelength Aq — 1334.5323 A 
which should provide, in conjunction with 7V(C IV), some indi- 
cation of the ionisation balance of C. For two of the C IV ab- 
sorbers in the present sample, at Zabs = 5.7238 and 5.8288 in 
line to SDSSJ103027.01-F052455.0 (Table |3j, the corresponding 
C II A1334.5323 lines are covered by the ESI spectrum of Pettini 
et al. (2003). In neither case is ClI absorption detected; the corre- 
sponding upper limits are NiC ll)/N(C IV) < 0.2 at Zabs = 5.7238 
and iV(C Il)/iV(C IV) ^ 0.05 at Zabs = 5.8288. Clearly these are 
high-ionisation absorbers for which it seems unlikely that C IV is 
only a trace ion stage of carbon. In future, it will be possible to 
quantify better the ionisation state of the metal-bearing IGM at 
these redshifts with observations targeted at other elements, such 
as Si which may be observable in three ion stages: Sill, Sim, and 
Si IV. 

Returning to Zigm, another reason why the value in eq. l |15t 
may underestimate the cosmic mean metallicity is that it does not 
include carbon in stars. The ensuing upward correction however, 
is unlikely to be more than a factor of ~ 2, unless current esti- 
mates (e.g. Bouwens et al. 2008) still miss a significant fraction of 
the cosmic star formation activity prior to z — 5.8 as argued, for 
example, by Faucher-Giguere et al. (2008). It is also intriguing to 
note that Schaerer (2002) proposed an efficiency factor as high as 
rj = 0.065 (nearly five times higher than the value used here in 
eq. [T6l > for metal-free stars. Additionally, the number of LyC pho- 
tons in the IGM may be boosted by a population of faint QSOs yet 
to be detected. 

Possibly, all of these factors contribute and, given all the un- 
knowns, one may argue that this kind of accounting exercise is pre- 
mature. Future work, aimed at improving the statistics of C IV ab- 
sorption at high redshifts and at identifying the galaxies from which 
these metals may have originated, will undoubtedly place the above 
speculations on much firmer observational ground. 



8 CONCLUSIONS 

We have conducted the largest survey to date for intergalactic met- 
als at redshifts 5.2 < z < 6.2. Our sample consists of near-IR 
spectra of nine high-redshift QSOs covering a total absorption dis- 
tance AX = 25.1, seven times greater than the pilot observations 
by Ryan- Weber et al. (2006). Our main results can be summarised 
as follows: 

1. We detect three definite and one possible CiV doublets; 
the three definite cases all correspond to column densities of triply 
ionised carbon in excess of 10^'* cm^^. We detect two further C IV 
doublets in a tenth QSO which is not included in the final sample 
as it shows evidence for mass ejection (a BAL QSO); in any case 
these two systems are probably associated with the QSO environ- 
ment rather than being truly intergalactic. 

2. From the three definite detections, we deduce a comoving 
mass density ficv = (4.4 ± 2.6) x 10~^ at a mean (z) = 5.76. 
This value represents a drop in Q,c iv by a factor of 3.5 compared 
to z ^ 4.7 at the 4(j significance level, after the different levels of 
completeness of optical and near-IR surveys are taken into account. 

3. If the column density distribution of CiV absorbers re- 
mained invariant at z > 4.7 in its shape and normalisation as de- 
termined at 1.5 ^ z ^ 4.5, we would have expected to find lOlg 
C IV doublets in our data, instead of only three, adding to the ev- 
idence for a decrease in the frequency of intergalactic C IV as we 
approach z ^ 6. 

4. Thus our data suggest a relatively rapid build-up of inter- 
galactic Civ over a period of only ~ 300 Myr, from ~ 0.9 to 
~ 1.2 Gyr after the Big Bang. Such a build-up could reflect the 
earlier rise of cosmic star formation activity from z ^ 9 indicated 
by galaxy counts and/or an increasing level of ionisation in the 
IGM. Numerical simulations of galaxy formation which include 
momentum-driven galactic outflows are consistent with our find- 
ings and suggest that both factors — a genuine build-up of metal- 
licity and higher ionisation — contribute to the rise in f^civ from 
z = 5.8 to z = 4.7. 

5. If the value of flc n we have derived is representative of 
the IGM at large, it corresponds to a metallicity Zigm ^ (9 ± 
5) X 10^^ Zq. The accumulation of this mass of metals in the IGM 
in turn implies that only about one Lyman continuum photon per 
baryon had been emitted by early-type stars prior to z — 5.8. This 
background is clearly insufficient to keep the the Universe ionised, 
suggesting that more carbon is present at these early epochs, in 
undetected ionisation stages, low column density systems, stars, or 
a combination of all three. 

6. All of these estimates are of necessity still rather uncer- 
tain. However, with the forthcoming availability of X-shooter on 
the VLtCI as well as planned improvements in the performance of 
near-IR spectrographs on other large telescopes, it should be pos- 
sible in the near future to make full use of metal absorption lines 
in QSO spectra to trace early star formation and the consequent 
reionisation of the IGM. 
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